Steady-state spectrum of Kelvin waves on a quantized vortex at finite temperatures 
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We determine numerically the steady-state spectrum for Kelvin waves on a superfiuid 4 He vortex that is 
driven by shaking its end points and damped by mutual friction. The relaxation rate towards a steady state is 
determined by the mutual friction. We have reached the zero temperature limit where the steady-state spectrum 
and the vortex length become independent of the temperature. Our drive with pointlike pinning sites produces 
a characteristic spectrum with \w m \ °= n-C^ , where 7] ~ 1.88. The spatially sharply peaked drive and the fact 
that we are fixing the oscillation amplitude, not the power, causes that there exists no high k cutoff from mutual 
friction even when mutual friction is large. This spectrum is generic to the pinned boundary conditions used. 
Without pinning the spectrum is sensitive to the drive. The spectrum even depends on the amplitude of the initial 
spectrum if vortex is allowed to relax on its own. Therefore, in real systems the Kelvin spectrum is not unique 
and depends on external conditions. 



I. INTRODUCTION 

During the last few years a central issue in the research on 
quantum turbulence (QT) has been the free decay of a vortex 
tangle. An important part of the decay at low temperatures is 
expected to be due to the Kelvin-wave cascade that works on 
length scales smaller than the intervortex distance. Due to cir- 
culation quantization the classical Kolmogorov decay cannot 
work at these scales. 

It is expected that at the lowest temperatures, where the dis- 
sipation due to mutual friction (m.f.) is negligible, the Kelvin 
cascade pushes the energy down to length scales of order the 
vortex core where it can be dissipated by phonons. How- 
ever, the Kelvin-wave cascade works also at finite tempera- 
tures with nonzero mutual friction but now there should exist a 
cut-off introduced by the m.f. Here we consider nonzero tem- 
peratures, since by doing so we avoid the difficulty to model 
the dissipation at the smallest scales. 

Kelvin waves (KWs) are helical distortions on a vortex 1 . 
Coupling between different Kelvin modes results in a cascade 
from large to smaller scales. Currently there are two main 
competing theories that predict a slightly different exponent 
for the Kelvin-wave spectrum 2,3 . There are also numerical 
calculations 4 " 7 , but the correct spectrum is still under active 
debate 8 " 14 . In the latest paper by Sonin 14 both theories are 
being criticized but the correct exponent for the Kelvin spec- 
trum is argued to be given by Kozik's and Svistunov's model, 
provided that the Kelvin amplitude is small enough. 

Here we calculate the Kelvin spectrum in a experimentally 
feasible situation where the Kelvin waves are exited by phys- 
ically shaking the vortex. We are most interested in the zero 
temperature limit. Our results indicate that the spectrum de- 
pends only weakly on mutual friction. Therefore we expect 
that our results are also valid in the limit T — >• 0. 

We also illustrate, by changing the drive, that the obtained 
spectrum is very generic, provided that the vortex end points 
are kept pinned. Without pinning the spectrum is much more 
sensitive to the drive. Also, if the vortex is allowed to relax on 
its own, the spectrum at later times depends on the amplitude 
of the initial spectrum. 



II. METHOD 

To obtain and accurately determine the steady state en- 
ergy spectrum one should be able to simulate a large range 
of scales. Even nowadays this is quite difficult when we have 
a large vortex tangle with hundreds of vortices. To obtain a 
steady state, one should additionally prepare a drive, that typ- 
ically works on large scales, and dissipation that is strongest at 
smallest scales. Since we are interested in the Kelvin cascade, 
we simplify the situation and consider only a single vortex in a 
configuration similar to that in Refs. [4-6]. Instead of simply 
removing the smallest scales, we damp the smallest scales by 
including mutual friction. To drive the Kelvin waves we force 
the end points of the vortex to move sinusoidally perpendicu- 
lar to the direction of the unperturbed straight vortex, i.e. x = 
Xq sin(o)f), y = 0, for points z = 0,L Z . We still assume periodic 
boundary conditions (b.c.) along the z-direction. This corre- 
sponds to a periodic grid of pinning points moving coherently, 
which is somewhat reminiscent of an oscillating grid where a 
vortex extends over several grid holes. The pinning points are 
considered to be small. The frequency of oscillation, ©, is 
adjusted to be close to the resonance frequency of the Kelvin 
wave, which for wave vector k m is given by 
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co m = -^[ln(2/ k m a)-y}. 
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Here K = h/m,4 « 0.0997 mm 2 /s is the circulation quantum, 
a w 10~ 7 mm is the vortex core size, and y = 0.57721 ... is 
the Euler constant. The above equation is valid for small am- 
plitude Kelvin waves in the long wavelength limit (k m a <C 1). 
Due to periodic boundary conditions the wave vector is dis- 
crete and given by k m = 2nm/L z , where m € Z. 

We use the vortex filament model where the vortex is con- 
sidered to be thin 15 . By adopting the scale separation scheme 
developed by Kozik and Svistunov 5 we are able to use much 
better resolution than in the earlier work by Vinen et al. 4 One 
should note that even if the scale separation scheme reduces 
the amount of work done (on average) per one time iteration 
from N 2 to N\nN, where N is the number of points on the vor- 
tex, the number of operations for a fixed time window scales 
like Af 3 \nN (without scale separation the work is proportional 
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to N 4 ). This is because the numerical time step scales like 
(L/N) 2 , where L is the vortex length. The scale separation 
routine has been criticized since the long and short scale dy- 
namics are partly disconnected and therefore it could give in- 
correct spectrum 10 . We have verified that in our case the ob- 
tained spectra with and without the scale separation scheme 
are identical when using a "medium" resolution (A^=1024). 

Our numerical scheme is standard practice in the vortex fil- 
ament formulation except for the addition of the scale sepa- 
ration scheme. We use the standard 4th order Runge-Kutta 
method and adapt many of the discretization methods de- 
scribed in the Ph.D. Thesis of Aarts 16 . The superfiuid velocity 
at point s is given by the Biot-Savart equation: 

K t (si -s) X dS\ 

\K J |S1-S| 3 

Here we have removed the singularity at Si = s by dividing 
the integral into a local and non-local term as in Ref. [15]. 
Vectors s' and s" in the local term are the tangent and the prin- 
cipal normal at s and the derivation is with respect to the arc 
length 4 . Terms Z_ and l + are the lengths of two adjacent line 
elements connected to the point s after discretization. In the 
non-local term the integral is now over the remaining config- 
uration. The numerical factor e 1 / 2 only takes into account the 
velocity distribution inside the core and is not essential here. 
At finite temperatures a vortex point moves at a velocity that 
is given by 15 

v L = v s + as' x (v n - T s ) - a's' x [§' x (v n - v g )] , (3) 

where a and a' are the temperature dependent mutual friction 
coefficients. For simplicity, we set a' = which is appropri- 
ate for 4 He at low temperatures. We additionally assume that 
the normal fluid is at rest, v„ = 0. The length (or period) of 
the unperturbed vortex is taken to be L, = 1 mm. The oscil- 
lation amplitude of the vortex end point is kept small enough 
to avoid reconnections and to keep the configuration such that 
we can describe it in terms of coordinates x(z) and y(z) as a 
sum of different Kelvin modes: 

x(z)+iy(z) =^w m expQ2nmz/L z ) . (4) 

111 

During simulations the Fourier components w m are deter- 
mined by using a FFT Since the vortex filament model does 
not automatically keep the points equally distributed along the 
z-axis we unify the grid before Fourier transformation by us- 
ing a second order interpolation. The error here is negligible 
since the difference from the uniform grid is small, provided 
that the grid unification is done frequently enough. 

The power dissipated due to mutual friction can be esti- 
mated by calculating P mf = /(v L • fmfMI = - /(vL • fu)d^, 
where f m f is force on a vortex due to mutual friction and 
is the Magnus force, both given per unit length (see Refs. [15 
and 17]). Assuming that v„ = the power dissipated becomes: 

P mf = -ap s K / |s' x \ s \ 2 d£, . (5) 




FIG. 1. (Color online) Steady-state Kelvin spectrum, averaged over 
the oscillation period of the drive. The drive is adjusted to be near 
mode in = 10 (a> = (Om) with amplitude xq = 0.0005 mm. The red 
curves are for a = 0.1, blue for a = 0.01, and magenta (in many 
places overlapping with blue) for a = 0.001. The simulations have 
been done using 256 (dash-dotted), 1024 (dashed), and 4096 (solid) 
points, which correspond to different cutoffs. The black dashed lines 
correspond to fits for a = 0. 1 and a = 0.01. The resulting slopes for 
\w m \ x m ^ are I? = 1-88 and 1.86, respectively. Inset: Steady-state 
vortex length L sat , also averaged over the oscillation period of the 
drive, as a function of the mutual friction parameter a. 



HI. STEADY STATE SPECTRUM 

Our simulations are started from an initial state with an 
unperturbed vortex without Kelvin waves. The drive is then 
switched on and simulations are continued until a steady state 
is achieved. The results are summarized in Fig. 1 where the 
steady-state Kelvin spectrum is plotted using different resolu- 
tions and different values of mutual friction. The spectrum is 
symmetric with respect to modes m and — m, therefore only 
positive modes are plotted. At least our "low" and "medium" 
resolution calculations indicate that the absolute value for the 
Kelvin amplitude \w m \ depends linearly on the drive ampli- 
tude. This was verified by increasing and decreasing the drive 
amplitude by a factor of five. In this respect our results differ 
from Ref. [4], where the Kelvin amplitude was observed to be 
independent of the drive amplitude. In our case dissipation is 
only due to mutual friction. Numerical dissipation (or noise) 
is quite small. For example, the kinetic energy is well con- 
served at T = 0. The largest error comes from the calculation 
of the m.f. dissipation at scales that are close to the numerical 
resolution limit m res , where the Kelvin-wave dispersion rela- 
tion is no longer accurately modeled, as noted in Ref. [5]. A 
simple estimation gives that co„, 1 — cos(7rm/ra res ), which 
becomes flat when m ~ m res and increases the timescale for 
modes near the resolution limit by a factor of 7T 2 /4 when 
compared with Eq. (I) 18 . Naturally, all the dissipation com- 
ing from scales smaller than our resolution is also omitted in 
the calculations. Therefore the calculation of P m f underesti- 
mates the real input power. At zero temperature in Ref. [4] 
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FIG. 2. (Color online) Transient response in vortex length L(t) as a 
function of time when the drive is switched on. The drive amplitude 
is xq = 0.0005 mm and its frequency is the resonance frequency of 
mode m = 10. The high frequency oscillations correspond to the 
frequency of the drive. The relaxation rate towards a steady state 
is mainly determined by mutual friction. The inset illustrates the 
additional oscillations, which become visible when a = 0.001. 
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the situation is more complicated in this respect since the de- 
termination of the input power is difficult. The amplitude of 
the oscillating superfluid velocity (which was used in Ref. [4] 
to drive the KWs) does not tell how much energy is really 
injected. In the pure Kelvin-wave cascade the absolute ampli- 
tude of the spectrum should still depend weakly on the amount 
of energy cascaded, e. Kozik-Svistunov's theory states that 
n m = \w m \ 2 °< e 1 / 5 , while L'vov-Nazarenko's formalism pre- 
dicts that «„, cx: g 1 / 3 . 

As noted above, the timescales for modes near the resolu- 
tion limit are longer than what Eq. (1) states. One might ex- 
pect that this explains why the spectrum bends up near m res . 
However, this is not the case since pinning three neighboring 
points (forcing these three points to move coherently), change 
the spectrum near m KS , while no effect is seen at lower m val- 
ues (Fig. 4). Therefore the spectrum near the numerical res- 
olution limit depends on the detailed structure of the pinning 
site. To consider length scales of order (and smaller) than the 
pinning site, one should take into account the geometry more 
precisely. 

The obtained spectrum with 77 w 1 .88 is close but not equal 
to the theoretical prediction by L'vov and Nazarenko 3 , where 
77 = 1 1/6 w 1 .833. This could indicate that the energy cascade 
is non-local and dominated by different-scale interactions. It 
is more likely that the approximate match of the exponent is 
accidental and specific to our case with pinning sites. This 
is discussed in more detail in Sec. V, where we analyze the 
effect of the drive and the boundary conditions. 



FIG. 3. (Color online) Time development of the Kelvin amplitudes 
after switch-on of the drive. The drive amplitude is xq = 0.0005 mm 
and its frequency is near the resonance frequency of mode m = 10. 
Left column of panels is for modes 10, 11, and 12 and right column 
for modes 20, 30, and 50. Note the different time intervals and scales. 



IV. TRANSIENTS AND RELAXATION TOWARDS THE 
STEADY STATE 

Due to the nature of the drive, the steady state is not a 
state where the vortex length and the amplitudes of the Kelvin 
waves are constant, but a state where they oscillate around 
their average values at a frequency determined by the drive 
frequency. These average values for the Kelvin amplitudes 
are plotted in Fig. 1. The oscillations are illustrated in Figs. 2 
and 3. The lower the temperature (smaller mutual friction pa- 
rameter a), the smaller is the oscillation amplitude (relative 
to the average value). Unfortunately, a smaller a also means 
a longer transient time to reach the steady state. The relax- 
ation time towards a steady state is inversely proportional to 
the mutual friction parameter a. Therefore we are mostly lim- 
ited to values a > 0.001. However, the simulations also indi- 
cate that the spectrum and the vortex length saturate below 
a = 0.01, which we interpret as the "zero temperature limit". 
Naturally, the value of a, below which the spectrum does not 
change any more, should depend on the scales considered. If 
the drive frequency would be much higher (corresponding to 
shorter Kelvin waves), one could expect that a smaller value 
of a would be required to reach the "zero temperature limit". 

During the relaxation towards the steady state one may ob- 
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FIG. 4. (Color online) Kelvin spectrum with pinned end points when 
N = 1024 and L z = 1 mm. The lower (blue) solid curve is the steady 
state spectrum for a = 0.01 and xq = 0.0005 mm, presented already 
in Fig. 1. Similar lower (blue) dash-dotted curve is calculated using 
the same parameters, but driving not only one, but three neighbor- 
ing points with the same amplitude. The upper (green) dash-dotted 
curve is for a vortex that is driven by oscillating superfluid velocity of 
the form v s = vocos(fcioz— (0\Qt)x with v'o = 0.05 mm/s, also using 
a = 0.01. The upper solid (red) curve describes the time averaged 
(over 0.25 s) spectrum after a vortex, with initial Kozik-Svistunov 
spectrum (rj = 1.7), has decayed 1.5 s with a = 0.001. The upper- 
most (black) dashed line illustrates the initial spectrum. The modes 
m = ±1 are modified to satisfy that x = y = at z = 0,L Z . 



serve oscillations at frequencies that are not related to the 
drive (Figs. 2 and 3). These oscillations are visible in the vor- 
tex length but also in the mode amplitudes when the mode is 
close to the drive mode. Oscillations with period T ~ 0.18 s 
are dominant for a = 0.01 and a = 0.001, but barely notice- 
able for a = 0.1. For a = 0.001 one may notice additional 
beating with a very similar frequency (inset of Fig. 2). It is 
not obvious whether these oscillations are related to Kelvin 
frequencies or not. The Kelvin period of mode m = 3 is 0.17 
s, but the neighboring modes are already quite far (periods re- 
lated to modes 2 and 4 are 0.37 s and 0.096 s, respectively). 
Generally, it seems that more oscillations, with different fre- 
quencies, become visible as a decreases and their amplitudes 
increase. It is possible that in the limit of zero temperature 
these oscillations grow such that a steady state becomes im- 
possible. 



V. IMPORTANCE OF THE BOUNDARY CONDITIONS 
AND THE EFFECT OF THE DRIVE 

The obtained spectrum is universal, independent of the 
drive, as long as the vortex end points are kept pinned. This is 
illustrated in Fig. 4, where the same spectrum appears when 
we drive the vortex by oscillating superfluid velocity of the 
form y s = t'ocos(£ m z— CO m t)x, similar to one used in Ref. [4] 
(while keeping the end points fixed at x = y = 0), or allow 



FIG. 5. (Color online) Kelvin spectrum without pinning, but still 
assuming periodicity, when N = 1024 and L- = 1 mm. The (black) 
dashed line is again the initial spectrum, with T] = 1.7, for the decay- 
ing vortex and the lower (red) solid curve describes the time aver- 
aged (over 0.25 s) spectrum after the vortex has decayed 1.5 s with a 
= 0.001. The upper (blue) solid line describes the steady state spec- 
trum when the vortex is driven in the &-space by periodically am- 
plifying the lowest modes |m| < 10 using a spectrum with r\ = 1.7 
and amplitude \w ±l \/L z = 0.035 with a = 0.002. The (blue) dash- 
dotted curves illustrate how a small increase in the drive amplitude 
(to |vv±i|/L z = 0.053) causes that the configuration "explodes" such 
that soon the Fourier presentation fails. 



the initial (e.g. Kozik-Svistunov) spectrum to relax towards a 
straight vortex. 

The situation is completely different if we omit pinning 
and allow the end points to move freely (within the limita- 
tion of the periodicity). In this case a high-fc mutual friction 
cut-off appears and practically all the modes above the drive 
frequencies are exponentially damped. Few example spectra 
are shown in Fig. 5. We were not able to observe a steady state 
spectrum of type \ w m \ mT^ . Decreasing the mutual friction 
damping or increasing the drive typically lifts up the high-£ 
part of the spectrum, but at some point even a small change 
in the drive or in the damping changes the spectrum from the 
exponentially damped one to a case where the Fourier pre- 
sentation fails. This is illustrated in Fig. 5, where we drive the 
vortex in the fc-space by amplifying the modes \m\ < 10 (with- 
out changing the phase) periodically such that these modes 
have the Kozik-Svistunov spectrum with, 77 = 1.7. More nu- 
merical points with much wider range of scales could perhaps 
help to avoid the over-amplification of the highest modes (by 
increasing dissipation), but since the required work for a fixed 
time window scales like N^lnN this becomes extremely time 
consuming to check. 

The strong difference in spectra with and without pinning 
sites can qualitatively be understood in the following way. 
Pinned end points cause a strong boundary for the Kelvin 
waves. KWs cannot travel freely but are typically reflected 
from the pinning site. This is similar to a case with plane 
boundary which prevents flow through it. Since we are con- 
sidering small pinning sites we are not forcing the vortex to 
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terminate perpendicularly to the pinning site but fix it loca- 
tion. This implies a small cusp in vortex, which means that 
more small scales structures appear. 

Next we will shortly describe results from simulation that 
mimic the calculations done by Kozik and Svistunov 5,6 . Here 
the mutual friction is omitted and the highest fc-modes are pe- 
riodically damped (exponentially) to mimic the energy sink. 
Pinning is omitted and periodic boundary conditions are used. 
The initial spectrum is L'vov-Nazarenko spectrum with 77 = 
11/6 and we follow the time development using different am- 
plitudes for the initial spectrum. Figure 6 summarizes these 
results. 

As one can see, the spectrum depends on the amplitude of 
the initial spectrum, as predicted by Sonin recently 14 . The 
larger the initial amplitude, the smaller 77 is observed. The 
smallest value for 77 = 1.41 is already so small that the 
vortex length would diverge without a cutoff near the core 
size. It is also a bit smaller than the lower limit predicted 
theoretically 14 . These values for 77, for a given amplitude, are 
only approximate since the spectrum still slowly changes with 
time. For example, the transition region from the low ^-region 
(having the initial spectrum) to the high fc-region with differ- 
ent 77 is quite well determined by equating the iteration time 
with the Kelvin period. The exponent 77 in the high ^-region 
still has a weak time dependence. To get rid off all the tran- 
sients one should continue iterations much further. However, 
without a drive the true steady state is difficult to reach. In 
any case, these simulations illustrate that the Kelvin spectrum 
is not unique, but depends on the amplitude. Theoretical pre- 
dictions by L'vov-Nazarenko and Kozik-Svistunov were done 
in the zero amplitude limit, which might be difficult to realize 
experimentally. 

The dash-dotted curves in Fig. 6, where the modes \m\ < 20 
are initially absent, also illustrate that spectrum is not simply 
determined by the amplitude but is also affected by the low- 
est modes. This would indicate that the cascade is non-local 
(and affected by the large scale curvature), supporting in this 
respect the theory by Lvov and Nazarenko. In any case, the 
spectrum is not as universal as one would expect from papers 
[5 and 6]. 



VI. ABSORPTION FROM MUTUAL FRICTION 

At high temperatures vortex motion is typically laminar and 
dissipation is only due to mutual friction. At low tempera- 
tures turbulence becomes more common and dissipation in- 
creases. In the zero temperature limit all the energy is cas- 
caded to scales of order the vortex core size and is dissipated 
by phonon emission ( 4 He). But at finite temperature, before 
the zero temperature limit is reached, the small scale Kelvin 
waves can greatly increase the m.f. dissipation when com- 
pared to the laminar case. This was seen in the simulations 
with oscillating pinning points (Sees. Ill and IV) where the 
dissipated power was strongly affected by the resolution used 
in the simulations. In order to see how fast the power dis- 
sipated increases with the number of Kelvin modes we inte- 
grated Eq. (5) by populating a vortex of length L z = 1 mm with 
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FIG. 6. (Color online) Kelvin spectrum during the decay at T = 
when N = 8192, L- = 1 mm, and the modes with \m\ > 1000 are 
exponentially damped to mimic the energy sink at large k. The solid 
curves are time averaged (over t = 0.015. . .0.020 s) and ^-smoothed 
spectra when the initial spectra (dashed lines) is \w m \ = A|77i|~ r ', with 
7J = 11/6. From bottom to top the initial amplitude, A/L z , is 0.0707, 

0. 106, 0.141, 0.177, and 0.233. For visibility, the different curves 
have been shifted vertically by multiplying them by 1/10, l/s/W, 

1, \/T0, and 10, respectively. The dash-dotted curves are for A/L z = 
0.233 (x 10) and 0.354 (x 10^/10) when initially the modes \m\ < 20 
are absent. For these runs the exponential damping of the highest 
modes is done less frequently since the energy flow is smaller and 
less damping is required. Inset: Compensated spectrum fc 17 / 10 ^! 
together with fits (black dashed lines), which estimate that rj = 1.72, 
1.60, 1.57, 1.48, 1.41, respectively. 



Kelvin waves, using up to 131072 discretization points on a 
vortex. Periodic b.c. are assumed. Note that here we do not 
pin (nor drive) the vortex end points. We just populate the vor- 
tex with a a spectrum \w„, \ 71 and include Kelvin waves 
up to some maximum mode /Vk.w (Nkw m res ), which corre- 
sponds to a cutoff scale Cutoff = L z /Nkw. The phase of each 
mode is set randomly and the amplitude for modes m = ± 1 
are |w±i|/L z = 0.01/V2- The positive and negative modes 
are chosen to have equal amplitude, i.e. \w m \ = \w- m \. This 
choice is not important, e.g. by choosing a random distribu- 
tion between negative and positive modes changes the results 
only by few percent. 

We observed no saturation in the dissipated power for 77 < 
2, at least for Kelvin spectra reaching modes up to 7Vkw ~ 
8000, which corresponds to a cutoff scale of order ^cutoff ~ 0. 1 
jj,m. For 77 = 2.5 the dissipated power might saturate as seen 
in Fig. 7, where the power dissipated is plotted as a func- 
tion of /View for a few different KW-spectra. A simple esti- 
mate, using the local induction approximation (LIA) where 
v. v = (k / An)\ri(R / a)s' x s", gives that the power converges, 
in the limit of /Vkw — > °°> only when 77 > 5/2. 

The dissipated power depends linearly on the m.f. param- 
eter a and, for small amplitudes, quadratically on the ampli- 
tude of the Kelvin waves. For a wide range of scales one can 
see from Fig. 7 that for 77 < 2 the power dissipated due to m.f. 
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FIG. 7. (Color online) Power dissipated by mutual friction as a func- 
tion of the largest Kelvin mode Afcw for different spectra \w m \ <* 
m^ 1 ! . The dashed lines show the best fits according to Eq. (6) when 
the four highest data points (blue circles) in each set are omitted from 
the fit. 

is approximately given by 

PM = aC^ oiV (6) 

where C depends at least on the amplitude of the KWs. We de- 
termined that j3 ps 1.38 for the Kozik-Svistunov 2 spectrum (77 
= 17/10) and j3 ps 1.14 for the L'vov-Nazarenko 3 spectrum (77 
= 1 1/6). The above equation implies that, if P m f is kept fixed 
by the drive and the spectrum is assumed to change only by a 
different (sharp) cutoff scale ^cutoff- an y change in a from a\ 
to 0£2 results in a corresponding change in Cutoff from ^cutoff, 1 
to ^cutoff.2 which are related as 

l^ = f^y /p . (7) 

<=cutoff,l \ a lJ 

Since the exponent l/j3 is close to 1, any change in a will 
cause a similar change in the cutoff scale. 

If one could obtain reasonable estimates for the cutoff scale 
and for the amplitudes of the lowest modes, one could try to 
estimate 77 by measuring the power dissipation. In practice 
this might be difficult. 

In summary, our calculations for power dissipation indicate 
that the mutual friction can dissipate lots of energy. This is es- 
pecially true for the case of 4 He where the core size is atomic 
and one can estimate that Afcw ~ P/ a ~ 10 6 , where £ ~ 0.1 
mm is the average vortex spacing. 

VII. DISCUSSION 

The calculated spectrum in Sec. Ill should correspond to an 
experimental situation where the pinning site is pointlike. An 
interesting point of these simulations is that there appears no 



m.f. dissipation dependent cutoff at high k-values, even if we 
use a relatively high value 0.1 for a. This is also caused by 
the spatially peaked drive and the fact that we are fixing the 
oscillation amplitude and not the power. By improving the 
numerical resolution, the power dissipated by mutual friction 
(which estimates the input power) increases quite rapidly with 
the number of Kelvin modes present, that is determined by our 
resolution. One should remember that the number of modes 
is limited by a physical cutoff determined by the vortex core 
size. At smaller scales, at least, the type of dissipation should 
change from m.f. dissipation to something else (e.g. phonons 
in case of 4 He). Unfortunately, the vortex filament formula- 
tion does not allow us to approach this limit. 

The Kelvin-wave cascade is expected to be a rather inef- 
ficient method for transferring energy to smaller scales and 
thus one might expect a bottleneck to appear in the energy 
transfer rate if the input power is too large 19 . Unfortunately, 
the ideal case to observe the bottleneck effect, low tempera- 
ture and high drive, is more demanding numerically and typ- 
ically results in a failure to present the vortex configuration 
as a Fourier series. Perhaps this failure is an indication of 
the breakdown of the cascade and the accumulation of Kelvin 
waves, owing to the appearance of the bottleneck. 

Numerical calculation of the energy flux through different 
length scales is a challenging problem since the flux due to 
Kelvin-wave cascade is extremely weak. Without a proper 
estimate for the energy flux it is difficult to quantify which 
part (or fraction) of the spectrum is generated by the Kelvin- 
wave cascade and what is the effect of the drive or due to 
pinning. As illustrated above, pinning tends to give a partic- 
ular spectrum, independent of the drive. Therefore, the spec- 
trum cannot be much affected by drive, but can be dominated 
by the effects due to pinning. By shaking the end points one 
forces an additional small U-shaped contour on the vortex. 
For example, a parabolic shape would result in a spectrum 
with \w m \ m~ 2 , when interpreted as a sum of Kelvin waves. 
However, since same spectrum appears for a decaying vortex, 
where this U-shaped contour is absent, the spectrum must be 
affected by a small cusp appearing at the pinning point and 
that causes a discontinuity in the derivative of the vortex con- 
figuration s(4). This discontinuity also results in a spectrum 
of type \w m \ °< nT 2 . However, sharp cusps in the vortex con- 
figuration are unavoidable at low temperatures with vanishing 
mutual friction. They often result from reconnections that are 
expected to drive the Kelvin cascade. 



VIII. CONCLUSIONS 

We have determined the Kelvin spectrum in a realistic, ex- 
perimentally feasible situation by driving the Kelvin waves by 
shaking the vortex end points and damping the high £- values 
by mutual friction. Therefore our results provide an improve- 
ment over previous work. When the vortex end points are kept 
pinned the obtained spectrum at scales larger than the pinning 
site is quite robust: independent of the drive and quite insen- 
sitive to temperature. Without pinning the spectrum is not 
unique. For example, the Kelvin spectrum for a vortex that 
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is allowed to relax on its own at T — produces a spectrum 
that depends on the absolute amplitude of the initial spectrum. 
Therefore, it is expected that in real systems the Kelvin spec- 
trum will depend on external conditions. 
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